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I. INTRODUCTION 



The B c meson is the ground state of the b-c system. It is a quarkonium state consisting 
of heavy quarks with different flavours, and it lies below the threshold of the BD mesons. 
It cannot annihilate into gluons, so it can only decay weakly. The B c system offers a new 
ideal place for studying the weak decay mechanism of heavy flavors and testing of the quark- 
: favour- mixing sector of the Standard Model. The discovery of the B c meson by the CDF 
l| has demonstrated the possibility for the experimental investigation of this system. 

B c mesons are too massive to produce at e + e~ colliders operating near T(45), such as the 
U B factories" at SLAC and KEK. They can be produced in significant numbers in hadron 
colliders. The Large Hadron Collider (LHC) is scheduled to run in 2008. Due to its high 
collision energy and high luminosity, the expected large production rates at the LHC will 
allow for precision measurements of B c properties. The 5-physics, including the physics 
potential of the B c system, could be fully exploited at the LHC. Measurements of heavy 
quarkonia at the LHC will be valuable for a deep understanding of the production and 
decay mechanisms of the corresponding quarkonia involved. 

The B c meson has very rich weak decay channels which can be divided into three classes: 
(1) the b quark decay {b — > c, u) with c quark as a spectator, (2) the c quark decay (c — > s, d) 
with b quark as a spectator, and (3) the annihilation channel. Clearly, we can study the two 
heavy flavors b and c simultaneously with B c meson. For the case of Class (1), the precision 
determination of the Cabibbo-Kobayashi-Maskawa (CKM) Q matrix elements of \V c b\ and 
\V u b\ are interesting, where determinations of \V u b\ based on inclusive and exclusive channels 
have an uncertainty of 8% and 17% at present, respectively [3j. For the case of Class (2), 
although the heavy B s (or B u , Bj) in the final state has a considerable effect in reducing 
the phase space for c quark decay, the CKM matrix elements \V c b\ <C |V CS | is in favor of the 
c quark decay greatly. For the case of Class (3), the B c weak annihilation decays are CKM 
enhanced by \V c b/V u b\ 2 ~ 10 2 as compared to B u annihilation decays. All the three classes 
listed above are interesting. 

The B c meson decay has been widely studied in the literature due to some of its 
outstanding features. The earlier nonleptonic decays of B c meson has been studied in 



outstanding ieatures. ine earner nonleptonic dec 
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19j. The theoretical status of the B c 



meson was reviewed in [19(]. In this paper, we will concentrate on the B c — ► J/ifiir, r/ c ir 
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decays in Class (1) with QCD factorization approach. Let us outline a few reasons below. 



The LHCb detector has a high trigger performance, good resolution for interaction 
vertex, excellent particle identification for charged particles. So all final-state parti- 
cles are detectable for the B c — > J/ipK, r] c ir decays. In practice, the most constructive 
information is to extract a signal of J ftp (or rj c ) from the cascade decays [20] . Espe- 
cially, compared with the semi-leptonic decays where the neutrino momentum is not 
detected directly, the LHC can provide a clean signature for B c — ► J/ipn, rj c 7c due to 
the narrow-peak of J/ip (or rj c ) and the good identification of the charged pion. It is 
estimated that one could expect around 5 x 10 10 B c events per year at LHC, so AT- 
LAS would be able to record about 5600 events of B c — > J/ip7i in one year [l9|. Better 
signal-to-background ratios and larger yields make B c — > J/iprc, i] c tt decay modes to 
be the most prospective channels for measurements. 

The final states of B c — > J/ip 1T , f] c ^ decays have large momentum compared with the 
Class (2) decays (e.g. B c — > B + X decays), so the final state interactions can be 
neglected among the energetic particles. One can effectively factorize the hard and 
soft amplitudes. Further more, the behavior of interaction vertices modified by hard 
gluon corrections at large recoil can be calculated reliably with perturbation theory. 

For the exploration of CP violation, the key processes are nonleptonic B decays. 
The theoretical challenge is the calculation of the hadronic transition matrix el- 
ements (f\Qi\B), where Qi is the local four quark operator. Recently, several 
attractive methods have been proposed to study the "nonfactorizable" effects in 



hadronic matrix elements, such as the QCD factorization (QCD 



tive QCD method [22|, |23j, |24j], soft collinear effective theory 
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) [21], perturba- 



261 ] . The nonlep- 



;onic two-bod y Pm , B s decays have been studied in detail, for example, in Refs 
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39fl . It is found that with appropriate 



parameters, most of the theoretical predictions are in agreement with the present ex- 
perimental data. With the accumulation of data, the Standard Model can be tested in 
more detail. The potential B c decay modes permit us to over-constrain the angles and 
sides of the unitarity triangle of the CKM matrix. The B c — > J/ipir, t] c tt decays are a\ 
dominant within the framework of Operator Product Expansion (OPE). They do not 
have pollution from penguins and annihilation diagrams. Moreover, the coefficient of 



at can be well determined compared with other coefficients. So these decays could be 
used for determination of the CKM matrix element \V c b\. 

In the limit of infinitely heavy quark masses, the HQET Lagrangians possess spin 
symmetry. The B c decays to J/ip, rj c involve the heavy spectator quark. The spectator 
enters the heavy hadron in the final states with no hard gluon rescattering, because 
spin symmetry works only at recoil momenta close to zero. The transition form factors 
should be nonperturbative hadronic quantities. It is expected that the quark-gluon sea 
is suppressed for B c mesons and that the calculation of factorization can be applied for 
non-leptonic B c decay. Thus the important parameter is the factor ai which depends 

i y, y, 0, y, s 0, y 



on the normalization scale 19j. The authors of 
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151, 



16 . 1 171 . 118H have much concern for hadron transition form factors whose dependence on 

n 

the momentum transfer is not significant within the physical phase space [191 ]. rather 
than the factor a±. Moreover, for the exploration of CP violation in B c decays, it 
needs to know the strong phases in addition to weak phases. While only the real value 
of the 
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actor ai_at t 
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re fixed point is used in the previous work 14, 

nn 



il, llTj, |18(] . In this paper, we will consider the radiative corrections 
to the hadron matrix element to compensate the renormalization scale dependence of 
Wilson coefficients. In addition, some information of strong phases can be obtained. 

This paper is organized as follows: In section [III we discuss the theoretical framework 
and compute the decay amplitudes for B c — > J/ipn, t] c tc with the master QCD formula. The 
section III 1 1 is devoted to the numerical results. Finally, we summarize in section IIVI 

II. THEORETICAL FRAMEWORK 
A. The effective Hamiltonian 



Using the operator product expansion and renormalization group equation, the low energy 
effective Hamiltonian for B c — > J/ipir, i] c ir decays can be written as 40] : 

Keff = + C 2 {jl)Q 2 } + H.c, (1) 

where V c bV* d is the CKM factors. The expressions of the local tree operators are 

Qi = {c a b a ) v _ A {dpUp) v _ A Q 2 = {c a bp) v _ A (dpu a ) v ^ A (2) 



Note that there is no local penguin operators, thus no interference between tree and penguin 
operators, and no CP asymmetry for these decays. The coupling parameters Ci(fi) are 
Wilson coefficients which have been evaluated to the next-to- leading order (NLO). Their 
numerical values in the naive dimensional regularization scheme (NDR) are listed in Table 

m 



B. Hadronic matrix elements within the QCDF framework 

For the weak decays of hadrons, the short- distance effects are well known and can be 
calculated with perturbation theory. However, the nonperturbative long-distance effects 
responsible for the hadronization from quarks to hadrons still remain obscure in several 
aspects. But to calculate the exclusive weak decays of the B c meson, one needs to evaluate 
the hadronic matrix elements, i.e., the weak current operator sandwiched between the initial 
state of the B c meson and the concerned hadronic final states, which is the most difficult 
theoretical work at present. Phenomenologically, these hadronic matrix elements are usually 
parameterized into the product of the decay constants and the transition form factors based 



on the naive factorization scheme (NF) 4l|. However, the main defect of the rough NF ap- 
proach is that the hadronic matrix elements cannot compensate the renormalization scheme- 
and scale- dependence of Wilson coefficients. In that sense the NF's results are unphysical. 
Moreover, for the exploration of CP violation, information of the strong phases is needed, 
but the values of both the hadronic matrix elements and the Wilson coefficients are real 
within NF. These indicate that "nonfactorizable" contributions from high order corrections 
to the hadronic matrix elements must be taken into account. 

A few years ago, Beneke, Buchalla, Neubert, and Sachrajda suggested a QCDF formula to 



compute the hadronic matrix elements [M\M 2 \Oi\B) in the 



reavy quark limit, combining the 



211 ] . At leading order in the heavy 



hard scattering approach with power counting in 1/rrib 
quark expansion, the hadronic matrix elements can be factorized into "non-factorizable" 
corrections dominated by hard gluon exchange and universal non-perturbative part param- 
eterized by the physical form factors and meson's light cone distribution amplitudes. This 
oromising method has been applied to exclusive two-body non-leptonic B u , Bd, B s decays 



33, 



34 



351 ] . It is found that with appropriate parameters, most of the QCDF's predictions 



are in agreement with the present experimental data. In this paper, we would like to apply 
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the QCDF approach to the case of B c — > J/i/jtt, t] c tt decays. 

In the heavy quark limit m& ^> Aqcd, up to power corrections of order of the 



master QCDF formula for B c — > r] c 7i (or J/ipir) can be written as 21 j 



(v c 7r\Oi\B c ) = F Bc ~* Vc ( l dz H(z)§ 7T (z) (3) 

J o 

Through the QCDF formula Eq.([3"l), the hadronic matrix elements can be separated into 
short- distance part [hard-scattering kernels H(z)] and long-distance part. At leading order 
of a s , H(z) = 1. There is no long-distance interaction between pion and B c -rj c (or B c -J/i/j) 
system, so the NF's picture is recovered. Nonperturbative effects are either suppressed by 
l/nib or parameterized in terms of the transition form factors F Bc ^ r,c (or F Bc ^ J ^) and the 
meson's light-cone distribution amplitudes &(z). 

C. Discussions on soft and collinear divergence cancellation at one-loop order 

When we calculate the vertex corrections in the leading power of Agnn/mh, not only 



ultraviolet divergence emerges but also infrared divergence, as shown in 42]. The authors 



of 21( gave an explicit cancellation of soft and collinear divergences in vertex corrections 
for B — > Dtt decays. They demonstrated that the soft divergences cancel out in Figfjja) 
and Figfjrjb), and in FigJT](c) and Fig{T](d), respectively; the collinear divergences cancel out 
in FigfjiTa) and FigJTJc), and in Figfjjb) and Fig{T](d), respectively. For B c — > J/ipn, t] c ti 
decays, the soft and collinear divergence cancellations are the same as that for B — > Dtt 



decays in 2l| . because they are all b — >• c transition. Here, we will consider the non-zero 



mass of c quark, and give an explicit calculation of the Feynman diagrams Figfjja)-(d) to 



show the cancellation of the infrared divergence, as shown in [42 1. 

The two valence quarks in pion have transverse momentum as well as longitudinal mo- 
mentum. The quark and antiquark momenta in pion meson is 

k q = zq + k±, k q = zq - k± (4) 

where q ~ 0(mb) is the momentum of the emitted n meson; z and z = 1 — z are the 
longitudinal momentum fraction of quark and anti-quark, respectively; k± ~ O(Aqcd) is 
the transverse momentum. But in B c — > J/ipir, t] c 7r decays, the pion is energetic. So 
the dominant contribution to the pion wave function comes from configurations where two 
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valence quarks are hard. The contribution from soft regions at the endpoint is believed to 
be suppressed in the heavy-quark limit. For our case, the quark is on-shell, and the effects 
proportional to k± is power suppressed. 

Now we shall take into account the a s corrections to the hadronic matrix element. After 
a straightforward calculation in the NDR scheme as in 42j, we obtain 

<<&>(.) = ^(*-\(du)v-A\0}(JM(cb)v-A\B;} 

( /m? v r(l-e) r l 2„ 1N l2 2lnt a „ „ T ./*a-l\ , 

x { " B) "IST b + 7 (hlt " " + ln 4 " + I3T " 4lni « " M-V") + 5 



t n hlt n 



(I - taf 



+ 



^ £ /._^£■n/ 1 , -U 1 *° ln *° 
£ 1-t 



+ (^)V)T ( i + e ; 



+ 1 



(5) 



(Q 2 } (6) = g^<7r-|(dTi)v_ A |0>(J/^|(a)v-A|S-> 



X 



{+B) 



' - + -(lnt b -l) + ln\ + , 



(4tt) 



*6 



41nt fe -2Li 2 (^— -) +6 



(^)Vmi+«) 



- + 4 

£ 1 



t b \nt b 



+ 11 



(6) 



a s C F 



(Q 2 ) (c) = ^^(7r-|(^) v _ A |0)(J/^|(c6) y ^| J B c -) 



x 



m c Vr(l- e ) 



(47T) 



1 1 .-, 21n£ r / 
- + - (lnt c - 1 + ln 2 t c + f - 41nt c - 2Li 2 



/i 2 \£,. ,r4 



't- - 1 



4 

- + 4 



1 - t r 



(Q2) 



t c \nt c 
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} 



M = ^^^1(d«)y-4|0)(^|(a)^|B c -) 



(7) 



X 



m 2 v r(l-e) r l 2., l2 21nf 
+ - (lnt d -l ) + ln 2 t d + : 



(47T) 



41nt d -2Li 2 (^-^) +5 



1 / t d \n t d 

r.vn - 



1 



/i 2 \ E/ . s ,=■„,.. x rl trflnt 



+ (^) £ (4 7 r)T(l + .) 



£ 1 - t d 



+ 1 



(8) 



where e — * + for ultraviolet divergences, e — > for infrared divergences, r c = m c /mb, and 

t a = ^(l-r 2 ) (9) 



* 6 = ^(l-r c 2 ) 
t c = z(l - l/r 2 c ) 
U = z(l - l/r 2 c ) 



-h/r 2 c 



(10) 

(11) 
(12) 
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In Eq.([7j), the term of infrared divergences can be written as 

- 2 ~-l 2 



-(lnt c - 1) + 



(£)'( 



1 



2 -t a 



1) 



— + ~(lni a - 1) - -(lnr 2 + i 2tt) + 



(13) 



We can treat the infrared divergences in Eq.(l8l) in the same way. 

- 2 ---l 2 



'771 

(^)'(^)' 



4 + + 



1 2., 
-a +- (lnt 6 



l)--(lnr 2 + z27r) + - 



(14) 



So it is easy to see that the infrared divergences are cancelled in the summation of all 
diagrams in Fig{T] before we integrate over the momentum fraction variable z. That is to 
say that the decay amplitude is infrared finite in the soft region when k± is neglected. This 
can be understood as a technical manifestation of colour transparency. 
Within the modified minimal subtraction (MS) scheme, we get 

a s C F 



Air N 

tnlllt 



{n-\(du) V - A \U){J/iP\(cb) v „ A \B-) 



+ 



l-*o 

t d \nt d 



l-U 

At c \nt c 

l-t c 



— r, 



+ 31n(-f)-9 



.1-4) 2 

1 r t d \nt d 



+ 



+ 



l -t 

1 



td 



where 



+2(lnt a - lnt b )lnr 2 - f(t a ) + f(t b ) + f(t c ) - f(t d )} 
l2 2t\nt nT . { t- In 

/(*) = lnt + T37- 2Ll2 (^) 



(15) 



(16) 



1 - * v t 

From above we found that the expressions of Eq. f)15j) and Eq. fll6l) are consistent with the 
result [see Eq.(80) — Eq.(85)] in the previous works 21]. 



With Eq. (jl5p and Eq. fll6p . we can compute the "non-factorizable" contribution of the one- 
gluon exchange vertex correction no matter the distribution amplitude <fi(z) is symmetric or 
not. This is very important in principle. For instance, when kaon is ejected from b quark 
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decay, the contribution from the asymmetric part of kaon distribution amplitude must be 
taken into account, although the contributions from the asymmetric part might be very 
small numerically. 

From the above, we can see that the "nonfactorizable" interactions connecting the B c — ► 
J/i/j (or i] c ) transition and the emitted pion are dominated by hard gluon exchange. So we 
may perform the k± integration over the pion wave function. We define [2l| 

d 2 k ± e ikd ' x . t* \ r A 
16tt 3 ^/2W c v ' ; L r5 M \ a( 3 



(7r(g) |4(^)^(0)|0) 



AN, 



dz (bJz) 



75 



a/3 



Assuming that the distribution amplitude (p n (z) is symmetric, i.e. (p n (z) 
Eq. ffT5l) can be simplified as follows: 

(Q2/( a ) + ( 6 ) + ( c ) + ( d ) 

= ^^-\(du)v-A\o)mmv- A \B:) 



(17) 
(z), then 



t a lnt a t d lnt 



t n 



-6 



/ t<a 



tnlntr. 



+ 



1-t, 



a 1 — 

+ -( // ■ xn + 



(18) 



Now we present the results of the coefficients a± for the case of B c — > J/ipir, i] c ir decays. 
Similar results hold for all Class (1) (ax-dominant) b — > egg' (where g and g' are light quarks) 
transition processes. 



„(NLQ) „ 
(NLO) 2 «S^F^(LO) 

The vertex corrections Vare given as follows [note: here <p(z) = <p(z) is assumed]: 
• for the m c ^ case 

2 2 

V = 31n(^)+31n(^)-18- fdz^H^z) 

where 

1 f t d \nt d 1 



(19) 



3t a lnt a 3t d lnt d / t a lnt a 1 \ 1 / f 



t d y 



i-t. 



(20) 



for the m,. = case 



32 



42 1 



2 — 

V = 3{21n(^) -6-i7r + Jdz<j)(z)^-^\nz} 



(22) 
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D. The renormalization scale dependence 



From the expressions of the coefficient a\ in Eg. ([19]) . apparently, the renormalization 
scale dependence of the hadronic matrix element is recovered, i.e. {Q 2 {n)) = g(n)(Q2)- We 
expect that the recovered scale dependence from can cancel the scale dependence of 
the Wilson coefficients Ci(fi), at least in part. 

With the renormalization group equations for the Wilson coefficients at leading order 
(LO) approximation 




An 



7. 



(0)T 



where the matrix of the LO anomalous dimensions 7^ is 

/ 




(23) 



7 



(0) 



V 



-6/iVc 
6 



6 ) 

Q/N c ) 



(24) 



It is easy to find that when the contributions from higher order of a s are neglected, we have 



dln/i 



a 1 



dln/i 



CM 



CM 



+ o ± C^ Cf J_ v 



4tt N r 



dln/i 







(25) 



Eq. (1251) just show that the coefficient a\ is no longer scale-dependent at the order of a s . Of 
course, a\ is still scale-dependent beyond the order of a s corrections. This can be seen from 
both Figj2] and the numbers in Table [B In principle, we have to include contributions of all 
higher order of a s corrections to cancel the remaining scale dependence of ai. 



E. Gauge dependence 

To demonstrate the gauge dependence, firstly we would like to write down the one-gluon 
exchange contribution to B c — > J Nm matrix element of the operator Q2. Following the 
expression of Eq.(55) — Eq.(57) in [2l|, we have 



where I is the momentum of the internal gluon propagator, and 

j> c - /+m c r r j) b + /+m b 

(Pc — 2 — m c (Pb + 2 — m 



A iV) = 7 A2_™ 2 r + r V ^A2__ 2 7 (27) 



Md.^ - r^/-/^r (28) 
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Here T = 7^(1—75) ; Pb and p c are the momentum of the b quark and the c quark, respectively. 
In fact, the terms in Eg. (126]) are calculated using Feynman gauge. For an arbitrary covariant 
gauge, the internal gluon propagator should be written as |43j 

7?(^ + ^) (29) 

where 77 = is the Feynman gauge, 77 = —1 is the Landau gauge. 

Assuming that the external states are all physical and can be approximated as on-shell 
quarks in the leading order of KQco/^bi it is easy to verify that for arbitrary covariant 
gauge, the terms proportional to 77 in Eq. (T2"6]) have no contribution, i.e. the expression of 
Eq. (126T) is gauge independent. For instance, 

/ / A<*m 1 y /+ m c r 1 ^ A+ / + m b 
a h 1 (0 " h \vc-lY-m{ + r (p 6 + Z)2-mg /h 

, 2 Pc -l-l 2 2p b -l + l 2 , , , 

= " I '^73p r + r ^7+^ = < 30 > 

Hence the gauge dependence cancels out when adding the two terms in A\. We can treat 
the A2 in the similar way and get the result, lal^A^ = 0. More precisely, we find that the 
expression of Eq. (l26|) is gauge independent. 



III. NUMERICAL RESULTS AND DISCUSSIONS 

The expressions of decay amplitudes for B c — ► J/ipTc, rj c Tr within the QCDF framework 
can be written as 

A{B-^J/^-) = v^G F K«/X^ J/ V/,(WP>i ( 31 ) 
A(B^ Vc n-) = -i^V^UF^iml-mlJa, (32) 

The branching ratios in B c meson rest frame can be written as: 

BR(B c ^X c5 tc) = ^^L\A(B c ^X cS n)\ 2 (33) 

where X c5 denotes J /if) (or i] c ) for B c — > J/ifm (or i] c tt) decays, and the momentum \p\ = 
{m 2 Bc — m 2 Xc _) / (2rriB c )- The lifetime and mass for B c meson are [3j: rriB c = 6.286 ± 0.005 
GeV, and t Bc = 0.46lg;J| ps. 

To perform numerical calculations, we specify the input parameters as follows. 

11 



Nonperturbative hadronic quantities, such as meson decay constants, transition form fac- 
tors, and meson distribution amplitudes, appear as inputs. In principle, information about 
these hadronic quantities can be obtained from experiments and/or estimated theoretically 
by non-perturbative method, such as lattice calculations, QCD sum rules, etc. In this paper, 
we shall follow the notation of 4J| on hadron transition form factors. As a good approxi- 
mation in the heavy quark limit, we shall take their values at the maximal recoil point for 
discussion. From the values collected in TableJITt we can see that there exists large uncer- 
tainty on theoretical predictions for F^°~ yric and Aq c ^ J ^. To estimate the branching ratios 
for B c — > J/ipn, i] c tt decays, we would like to take F^ c ^ Vc = Aq°~~* J = 0.6. We shall use 
the asymptotic form of the pion in the calculation, i.e. <j) n (z) = 6zz. 

The Wolfenstein parameterization for the CKM matrix elements is used due to its ad- 
vantage in explicitly expressing the hierarchy among the CKM elements in terms of powers 
of small parameter A. With this notation, up to the order of 0(X A ), the CKM elements 
involved are V c b = AX 2 , and V u d = 1 — A 2 /2. The independent parameters of A and A have 
been well determined to the accuracy of 1% and 0.1%, respectively. Their values are [3] 

A = 0.818lg;{j?£, A = 0.2272±0.0001 (34) 

Note that the relative weak phase difference is zero at this approximation, so no CP- violating 
asymmetry is expected. 

Other input parameters are 



m u = rrid = m s = 0, m c = 1.25±0.09 GeV 

m b = 4.20±0.07 GeV f n = 130.7±0.1±0.36 MeV 



= 3096.916±0.011 MeV Q m vc = 2980.4±1.2 MeV 

If not specified explicitly, we shall take their central values as the default input. 
Our numerical results are listed in Table [I] and IIII1 

From the numbers in Table [H we can see that part information of strong phases is obtained 
by considering gluon radiative corrections to vertexes. Compared with the real part, the 
imaginary part of eti is small for B c — * J/ipn, r] c 7c. The reason is that the "nonfactorizable" 
effects are a s or AQco/mb suppressed within QCDF. Although the imaginary part is small, 
it will be very important for the exploration of CP violation for decay modes receiving 



contributions from both tree and penguin topologies. In fact, the numerical results of 
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421 ] indicate that the imaginary parts of most Oj (i = 3, • • •, 10) are the same order as 



the real parts, i.e. the "nonfactorizable" effects cannot be neglected for penguin-dominated 
decay modes. 

It should be noted that both the imaginary part arising at the order of a s and the real 
part of ai depend on the renormalization scale. But the scale dependence of d\ has been 
reduced compared with that of NF, at least at the order of a s . This can be seen from the 
Eq. fT25l) and FigfSJ In principle, high order of a s corrections could be calculated order by 
order. Perhaps the remaining scale dependence of ai could be cancelled by the contributions 
from higher order of a s corrections. 

The numerical results of the branching ratios for B c — > J/ip7i, i] c 7[ decays are given in 
Table IIHI From the numbers in Table IIII[ we can see that these two branching ratios are 
close to each other, if we assume Aq c ^ J ^ = F^ c ^ Vc . In fact, this point can also be seen 
from their relationship 



B c ^J/i> |2/_ 2 ^2 



BK{B C -^J/^) K c j/ T« c 



m 



(35) 



Some nonperturbative effects could be cancelled in the ratio of branching fraction. Using 
experimental data, the relation between the form factor Aq c ^ J ^ and F^ c ^ Vc is expected 
to be obtained from the ratio Eq. fl35|) . Once the nonperturbative parameters Aq c ^ j ^ and 
Fq c ^ Vc are fixed, we might get some information on CKM element V c b from Eq. (!3~TT) and 
Eq.([32D. 

In addition, form Table III we also see that the charm quark mass effect on the branch- 
ing ratios is small, and that uncertainties form renormalization scale is also small. The 
uncertainties from the CKM elements which is proportional to 2a a/ A are less than 5% [see 
Eq. (!34j) ]. Thus the remaining uncertainties are the hadron parameters, such as transition 
form factors. Within QCDF approach, the transition form factor is the nonperturbative 
input parameter which comes from the long-distance contributions, so it must be computed 
nonperturbatively or determined experimentally. From the numbers in Table [TTJ, we can see 
that there are large differences in the numerical values of transition form factors Aq c ^ J ^ , 
F c for various theoretical approaches. So at the scale of \i = mi,, the branching ratio 
can be rewritten as 

Bn{B c ^J/i,7i) = 0.119%x ' ° Q 36 1 (36) 
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BTZ(B c ^ Vc 7r) = 0.128%x |F ° J ^ (37) 

0.36 

The effects of transition form factors on the branching ratios are displayed in FigfS] a , in- 
cluding the uncertainties from the CKM elements, quark masses and the renormalization 
scale. (Here the uncertainties from Aq C D/ m b corrections is not included, but it is assumed 
that its contribution is power suppressed and should not be large 32J) Clearly, the largest 
theoretical uncertainty comes from the transition form factors, which is of a nonperturba- 
tive nature within the QCD framework. Here we can not compute these form factors within 
QCD factorization scheme. We just use them as input parameters. 



IV. SUMMARY AND CONCLUSION 



In this paper, we calculated the hadronic matrix element for b — > c transition at one- 
loop level under NDR scheme and in the heavy quark limit. Then we apply the master 
QCD formula to the would-be well detectable decay channels in experiments B c — > J/ipK, 
i] c 7T. The "nonfactorizable" vertex corrections are computed. We find that in the heavy 
quark limit, the "nonfactorizable" contributions dominated by hard gluon exchange can 
compensate the renormalization scale-dependence of the Wilson coefficients and are infrared 
safe and gauge independent. Finally, we give the branching ratios for B c — > J j tpn, f] c 7i decays 
using the asymptotic distribution amplitude of the pion. We find that the large theoretical 
uncertainties come mainly from the non-perturbative transition form factors. 
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TABLE I: Wilson coefficients Cj and a\ in the NDR scheme. 



1* 


NLO 


LO 


ATTT 1 


QCDF 


(GeV) 


Ci 




Ci 




a\ 


ax 

(m c = 1.25 GeV) 


ax 

(m c = 0) 


2.0 


1.147 


-0.304 


1.191 


-0.374 


1.046 


1.084 + i 0.034 


1.079 + i 0.036 


3.0 


1.109 


-0.236 


1.148 


-0.303 


1.031 


1.070 + i 0.024 


1.066 + i 0.025 


4.0 


1.088 


-0.196 


1.124 


-0.261 


1.023 


1.061 + i 0.019 


1.058 + i 0.020 


5.0 


1.074 


-0.169 


1.108 


-0.232 


1.018 


1.054 + i 0.016 


1.052 + % 0.016 


6.0 


1.064 


-0.148 


1.096 


-0.210 


1.014 


1.049 + i 0.013 


1.047 + % 0.014 


7.0 


1.056 


-0.131 


1.087 


-0.192 


1.012 


1.045 + i 0.012 


1.043 + i 0.012 


8.0 


1.049 


-0.118 


1.079 


-0.178 


1.010 


1.041 + i 0.011 


1.040 + i 0.011 



TABLE II: Values of transition form factors 



Reference 




A B c ^J/4> 


[4] 

m 

[44] 
[45] 
[46] 
[47J 
[14] 
[48] 
[18] 


0.170 ~ 0.687 
0.20±0.02 
0.23±0.01 

0.66 

0.76 

0.87 

0.47 

0.61 


0.156 ~ 0.684 
0.26±0.07 a 
0.21±0.03 6 
0.595 c 
0.69 d 
0.27 e 
0.40 
0.57^ 
^7+ 01 

u - ' -0.02 
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TABLE III: Branching ratios for B c — ► J/tpn, r] c ir (in the unit of %). 







BTl{B c - 


+J/ipir) 




BK{B C - 


-*»7ct) 




NF 


QCDF 


QCDF 


NF 


QCDF 


QCDF 






(m c = 0) 


(m c = 1.25 GeV) 




(m c = 0) 


(m c = 1.25 GeV) 


m 6 /2 


0.116 


0.123 


0.124 


0.124 


0.132 


0.134 


m 6 


0.111 


0.119 


0.119 


0.119 


0.127 


0.128 


2mb 


0.108 


0.114 


0.115 


0.116 


0.123 


0.123 





(a) (b) (c) (d) 

FIG. 1: Vertex corrections to hard-scattering kernel for b — > c decay at the order of a s . The upward 
lines represent the valence quarks of the emitted n meson. 



Compared the definitions of transition form factor of Ref. 41] with those of Ref.[6|], we can obtain their 
relationships at the maximal recoil point, 



.4 



B e ^J/tl> 




2m , 



2m , 



(38) 



where the transition form factor A 



o 



is defined in 



4l|, Fq andFf are defined in 6] , and their values 



are F A = 2.5±0.3 GeV" 1 , F A = -0.03±0.01 GeV" 1 . 
b Using the relationship of Eq.® and the values of = 2.0±0.2 GeV -1 , Ff = -0.024±0.002 GeV" 1 



44] 



44] 



c Using the relationship of Eq.{3S]| and the values of F A = 5.9 GeV" 1 , F A = -0.Q7J GeV" 1 
d Compared the definitions of transition form factor of Ref. [4l| with those of Ref. [46( , we can obtain their 
relationships at the maximal recoil point, 



A 







2m, 



2m 



'M±A 2 



(39) 



where A\ and A 2 are defined in [46(, and their values are A\ — 0.68, A 2 =fl-Q6. 
Using the relationship of Eq. (|3T))) and the values of A\ = 0.75 . A 2 = 1.69. 



f Compared the definitions of transition form factor of Ref. 41 j with those of Ref. 48|, we can obtain their 
relationships at the maximal recoil point, 



A 



B c ^>J/ip 




'J/4: 



2m 



(Ao - A+) 



(40) 



//</> 



where Aq and A + are defined in 481, and their values are Aq = 1.64, A 2 = 0.54. 
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1 . 1 1 0.05 




3456789 10 3456789 10 

(a) (b) 
FIG. 2: Dependence of the coefficient a\ [vertical axes, Re(ai) in (a) and Im(oi) in (b)] on the 
renormalization scale fi [horizontal axes, in units of GeV], with asymptotic light-cone distribution 
amplitudes (p w (x) = 6xx. The solid and dashed lines denote a\ for the case of m c = 1.25 GeV and 
m c = within QCDF framework, respectively; the dotted line denotes a\ within NF framework. 
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A B c ->J/tp TiBc^ric 
A) ^0 

FIG. 3: Branching ratios for B c — » J/i/jit, ?7 c vr versus the transition form factors Aq c ~ >j ^ ) pRc-^Vc 
in the QCDF approach, considering the mass of c quark. The solid lines are calculated with central 
values of inputs, while the bands denote the uncertainties from CKM elements and quark masses. 
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